use "C:\Users\andre\Desktop\Artigos\Artigo_Ens_Sup_2\Bases\Agregadas_Stata\Agregada.dta", clear

estimates drop _all


********************************Modelo para conclusão do EM

*Sem controles
logistic ens_medio b(first).ano##b(1).quintis_renda if (rec_mod==1) & ((ano==1992 | ano==2002 | ano==2022) | (ano==2012 & versao==1)) [pweight = peso] 

*Com controles

*1992
quietly logistic ens_medio b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==1992) & (rec_mod==1) & (((ano>=1992 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_renda_1992

quietly margins quintis_renda, post
estimates store mg_em_renda_1992

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_ens_medio_renda_1992

*1993
quietly logistic ens_medio b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==1993) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_renda_1993

quietly margins quintis_renda, post
estimates store mg_em_renda_1993

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_ens_medio_renda_1993

*1995
quietly logistic ens_medio b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==1995) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_renda_1995

quietly margins quintis_renda, post
estimates store mg_em_renda_1995

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_ens_medio_renda_1995

*1996
quietly logistic ens_medio b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==1996) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_renda_1996

quietly margins quintis_renda, post
estimates store mg_em_renda_1996

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_ens_medio_renda_1996

*1997
quietly logistic ens_medio b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==1997) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_renda_1997

quietly margins quintis_renda, post
estimates store mg_em_renda_1997

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_ens_medio_renda_1997

*1998
quietly logistic ens_medio b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==1998) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_renda_1998

quietly margins quintis_renda, post
estimates store mg_em_renda_1998

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_ens_medio_renda_1998

*1999
quietly logistic ens_medio b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==1999) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_renda_1999

quietly margins quintis_renda, post
estimates store mg_em_renda_1999

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_ens_medio_renda_1999

*2001
quietly logistic ens_medio b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2001) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_renda_2001

quietly margins quintis_renda, post
estimates store mg_em_renda_2001

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_ens_medio_renda_2001

*2002
quietly logistic ens_medio b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2002) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_renda_2002

quietly margins quintis_renda, post
estimates store mg_em_renda_2002

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_ens_medio_renda_2002

*2003
quietly logistic ens_medio b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2003) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_renda_2003

quietly margins quintis_renda, post
estimates store mg_em_renda_2003

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_ens_medio_renda_2003

*2004
quietly logistic ens_medio b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2004) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_renda_2004

quietly margins quintis_renda, post
estimates store mg_em_renda_2004

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_ens_medio_renda_2004

*2005
quietly logistic ens_medio b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2005) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_renda_2005

quietly margins quintis_renda, post
estimates store mg_em_renda_2005

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_ens_medio_renda_2005

*2006
quietly logistic ens_medio b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2006) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_renda_2006

quietly margins quintis_renda, post
estimates store mg_em_renda_2006

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_ens_medio_renda_2006

*2007
quietly logistic ens_medio b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2007) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_renda_2007

quietly margins quintis_renda, post
estimates store mg_em_renda_2007

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_ens_medio_renda_2007

*2008
quietly logistic ens_medio b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2008) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_renda_2008

quietly margins quintis_renda, post
estimates store mg_em_renda_2008

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_ens_medio_renda_2008

*2009
quietly logistic ens_medio b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2009) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_renda_2009

quietly margins quintis_renda, post
estimates store mg_em_renda_2009

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_ens_medio_renda_2009

*2011
quietly logistic ens_medio b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2011) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_renda_2011

quietly margins quintis_renda, post
estimates store mg_em_renda_2011

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_ens_medio_renda_2011

*2012
quietly logistic ens_medio b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2012) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_renda_2012

quietly margins quintis_renda, post
estimates store mg_em_renda_2012

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_ens_medio_renda_2012

*2013
quietly logistic ens_medio b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2013) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_renda_2013

quietly margins quintis_renda, post
estimates store mg_em_renda_2013

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_ens_medio_renda_2013

*2014
quietly logistic ens_medio b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2014) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_renda_2014

quietly margins quintis_renda, post
estimates store mg_em_renda_2014

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_ens_medio_renda_2014

*2015
quietly logistic ens_medio b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2015) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_renda_2015

quietly margins quintis_renda, post
estimates store mg_em_renda_2015

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_ens_medio_renda_2015

*2016
quietly logistic ens_medio b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2016) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_renda_2016

quietly margins quintis_renda, post
estimates store mg_em_renda_2016

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_ens_medio_renda_2016

*2017
quietly logistic ens_medio b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2017) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_renda_2017

quietly margins quintis_renda, post
estimates store mg_em_renda_2017

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_ens_medio_renda_2017

*2018
quietly logistic ens_medio b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2018) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_renda_2018

quietly margins quintis_renda, post
estimates store mg_em_renda_2018

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_ens_medio_renda_2018

*2019
quietly logistic ens_medio b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2019) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_renda_2019

quietly margins quintis_renda, post
estimates store mg_em_renda_2019

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_ens_medio_renda_2019

*2020
quietly logistic ens_medio b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2020) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_renda_2020

quietly margins quintis_renda, post
estimates store mg_em_renda_2020

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_ens_medio_renda_2020

*2021
quietly logistic ens_medio b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2021) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_renda_2021

quietly margins quintis_renda, post
estimates store mg_em_renda_2021

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_ens_medio_renda_2021

*2022
quietly logistic ens_medio b(1).quintis_renda b(first).sexo b(first).cor c.idade b(first).regiao b(first).local if (ano==2022) & (rec_mod==1) & (((ano>=1993 & ano<=2015) & (versao==0)) | ((ano>=2016 & ano<=2022) & (versao==1))) [pweight = peso] 
estimates store mod_ens_medio_renda_2022

quietly margins quintis_renda, post
estimates store mg_em_renda_2022

quietly nlcom (renda2_EM: _b[2.quintis_renda] / _b[1.quintis_renda]) (renda3_EM: _b[3.quintis_renda] / _b[1.quintis_renda]) (renda4_EM: _b[4.quintis_renda] / _b[1.quintis_renda]) (renda5_EM: _b[5.quintis_renda] / _b[1.quintis_renda]), post
estimates store nlcom_ens_medio_renda_2022


****************Tabelas

*OR
etable, estimates(mod_ens_medio_renda_1992 mod_ens_medio_renda_1993 mod_ens_medio_renda_1995 mod_ens_medio_renda_1996 mod_ens_medio_renda_1997 mod_ens_medio_renda_1998 mod_ens_medio_renda_1999 mod_ens_medio_renda_2001 mod_ens_medio_renda_2002 mod_ens_medio_renda_2003 mod_ens_medio_renda_2004 mod_ens_medio_renda_2005 mod_ens_medio_renda_2006 mod_ens_medio_renda_2007 mod_ens_medio_renda_2008 mod_ens_medio_renda_2009 mod_ens_medio_renda_2011 mod_ens_medio_renda_2012 mod_ens_medio_renda_2013 mod_ens_medio_renda_2014 mod_ens_medio_renda_2015 mod_ens_medio_renda_2016 mod_ens_medio_renda_2017 mod_ens_medio_renda_2018 mod_ens_medio_renda_2019 mod_ens_medio_renda_2020 mod_ens_medio_renda_2021 mod_ens_medio_renda_2022) cstat(_r_b, nformat(%4.3f)) stars( , attach(_r_b) ) mstat(N, nformat(%9.0g))

collect layout (coleq#colname#result[_r_b] result[N] result[aic] result[bic] result[chi2] result[df_m] result[p] result[ll] result[r2_p]) (etable_depvar#stars) (), name(ETable)
collect export "C:\Users\andre\Desktop\Artigos\Artigo_Ens_Sup_2\Resultados\Modelos.xlsx", as(xlsx) sheet(Ens_Med_renda_OR, replace) cell(A1) modify


*RR
etable, estimates(nlcom_ens_medio_renda_1992 nlcom_ens_medio_renda_1993 nlcom_ens_medio_renda_1995 nlcom_ens_medio_renda_1996 nlcom_ens_medio_renda_1997 nlcom_ens_medio_renda_1998 nlcom_ens_medio_renda_1999 nlcom_ens_medio_renda_2001 nlcom_ens_medio_renda_2002 nlcom_ens_medio_renda_2003 nlcom_ens_medio_renda_2004 nlcom_ens_medio_renda_2005 nlcom_ens_medio_renda_2006 nlcom_ens_medio_renda_2007 nlcom_ens_medio_renda_2008 nlcom_ens_medio_renda_2009 nlcom_ens_medio_renda_2011 nlcom_ens_medio_renda_2012 nlcom_ens_medio_renda_2013 nlcom_ens_medio_renda_2014 nlcom_ens_medio_renda_2015 nlcom_ens_medio_renda_2016 nlcom_ens_medio_renda_2017 nlcom_ens_medio_renda_2018 nlcom_ens_medio_renda_2019 nlcom_ens_medio_renda_2020 nlcom_ens_medio_renda_2021 nlcom_ens_medio_renda_2022) cstat(_r_b, nformat(%4.3f)) cstat(_r_lb, nformat(%4.3f)) cstat(_r_ub, nformat(%4.3f)) export("C:\Users\andre\Desktop\Artigos\Artigo_Ens_Sup_2\Resultados\Modelos.xlsx", as(xlsx) sheet(Ens_Medio_renda_RR, replace) cell(A1) modify)
